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Abstract 



We consider the decay of a heavy flavour into an inclusive hadronic state 
\q \ X of invariant mass mx small with respect to its energy Ex, mx Ex- 

The electron spectrum and the hadronic mass distribution in semileptonic 
b — > u decays, or the photon spectrum in b — > 37 decays, all require, close 
to their endpoints, a control over this region. This region is affected both 
by non-perturbative phenomena related to the Fermi motion of the heavy 
quark and by perturbative soft gluon radiation in the final state (Sudakov 
form factor). Fermi motion can be described by the shape function /(m*), 
Ph! which represents the distribution of the effective mass m* of the heavy quark 

at disintegration time. We perform a factorization with a simple technique 
in order to consistently separate perturbative from non-perturbative effects. 
We find that the shape function, contrary to naive expectations, is not a 
physical distribution, as it is affected by substantial regularization scheme 
effects, controlling even the leading, double-logarithmic term. It factorizes, 
however, the bulk of non-perturbative effects in lattice-like regularizations. 
Some non-perturbative effects are present in the coefficient function even at 
leading twist, but they are expected to be suppressed on physical grounds. 
Finally, we clarify a controversial factor of 2 in the evolution kernel of the 
shape function. 
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1 Introduction 



Nowadays there are many facilities that allow an accurate experimental study 
of heavy flavour decays. It is therefore becoming more and more important 
to improve the accuracy and the reliability of the theoretical calculations. In 
this paper, we study the properties of the decays of heavy flavour hadrons 
into inclusive hadronic states X "tight" in mass, i.e. with an invariant mass 
mx small with respect to the energy Ex'- 

mx < E x . (1) 
More specifically, we consider the situation where 



O {E x Aqcd), (2) 



O « 1 {E x » Aqcd) • (3) 



so that 

™* ~ o ( 
E 2 X V E x 

A formal definition of kinematics @ is the limit, in the heavy quark rest 
frame: 

E x -> oo, 
m x — > oo 



with 



x 



const. (4) 



E x 

The divergence of m x - even though slower than the one of E\ - implies 
that the final hadronic state can be replaced with a partonic one, i.e. that 
the use of perturbation theory is fully justified. Heavy flavour decays are 
characterized by three mass or energy scales: the mass of the heavy flavour 
ttlq, the energy Ex, and the invariant mass mx of the final hadronic state. 
Limit (f|) implies also the limit of infinite mass for the heavy flavour Q: 

m Q -> oo, (5) 

since uiq > Ex- Another consequence of is that 

T# - 0, (6) 
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i.e. we are in the so-called threshold region^. 

The study of these processes has both a theoretical and a phenomenologi- 
cal interest. On the theoretical side, in heavy decays the infrared perturbative 
structure of gauge theories - the Sudakov form factor |], Q - enters in a rather 
"pure" form, owing to the absence of initial state mass singularities. On the 
phenomenological side, the computation of many relevant distributions re- 
quires a good theoretical control over the region ([!]). As examples, let us 
quote the electron spectrum dT/dE e close to the endpoint E e < [|3| 
and the hadron mass distribution dT/dmx at small mx H in semileptonic 
b — > u decays, such as 

B -> X u + e + v, (7) 

or the photon energy distribution dT/dE^ close to the endpoint E 1 < mg/2 
in b — > S7 decays. For the electron or photon spectrum, the region (Ij) is 
involved because the requirement of a large energy of the lepton or of the 
photon pushes down to zero the mass of the recoiling hadronic system. As 
is well known, the above mentioned distributions in (|7]) allow an inclusive 
determination of the CKM- matrix element \V u b\ ||, while a large photon 
energy in the rare decay, E n > 2.1 GeV is required to cut experimental 
backgrounds. 

In general, the dynamics in region (JI]) is rather intricate as it involves 
an interplay of non-perturbative and perturbative contributions. These are 
related to the Fermi motion of the heavy quark inside the hadron and to the 
Sudakov suppression in the threshold region (^), respectively. Even though 
these two effects are physically distinguishable and are treated as indepen- 
dent in various models ||, they are ultimately both described by the same 
quantum field theory, QCD. Therefore the problem arises of describing them 
consistently, i.e. without double countings, inconsistencies, etc. Our idea is 
to subtract from the hadronic tensor encoding all QCD dynamics, 

W, v = J2(HQ\J+\X){X\J^\H Q )5 4 {p B -q-px), (8) 
x 

each of the perturbative components - associated with the Sudakov form 
factor and with other short- distance corrections - to end up with an explicit 

1 The converse is not true: limit (Q) does not imply limit (|J) (we thank G. Veneziano 
for pointing this out to us). 
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representation of the non-perturbative component. In eq. (|8|), we have defined 



J M (x) = q(x)Y^Q(x), (9) 

where V ^ is a matrix in Dirac algebra^, px is the momentum of the final 
hadronic jet, and Hq is a hadron containing the heavy quark Q. The non- 
perturbative component is identified with an ultraviolet (UV) regularized ex- 
pression for the structure function, or shape function, in the effective theory. 
The shape function f(k + ) has been introduced using the Operator Product 
Expansion (OPE) and can be defined as || 

f{k+) = (H Q I h\ 5(k + - iD + ) h v | Hq), (10) 

where h v is a field in the Heavy Quark Effective Theory (HQET) with 4- 
velocity v; D + is the plus component of the covariant derivative, i.e. D + = 
D° + D 3 . The shape function represents the probability that the heavy quark 
has a momentum tubV + k' with a given plus component k' + = k + . This 
function can also be interpreted (see section [Oj) as the probability that Q 
has an effective mass 

m* = m B + k + (11) 

at disintegration time. The renormalization properties of the shape function 
have also been analysed [7|, R BL [111]. Because of UV divergences affecting 
its matrix elements, f(k + ) needs to be renormalized and it consequently 
acquires a dependence on the renormalization point fi: f(k + ;fi). The non- 
perturbative information about Fermi motion enters in this framework as the 
initial value for the /i-evolution. The shape function can be extracted from 
a reference process and used to predict other processes, analogously to the 
parton distribution functions in usual hard processes such as Deep Inelastic 
Scattering (DIS) or Drell-Yan JTT]]. In principle, it can also be computed 



with a non-perturbative technique, for example lattice QCD [|T2| 



Our approach aims at a deeper understanding of perturbative and non- 
perturbative effects with respect to the standard OPE in dimensional regu- 
larization (DR). We compare different regularization schemes and find that 
the factorization procedure is substantially scheme- dependent. By that, we 
mean a much stronger scheme dependence than the usual one, typically re- 
lated to the finite part of one-loop amplitudes, corresponding in DR to a 



2 For the left-handed currents of the Standard Model, = 7 M (1 — 75) 
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replacement of the form 1/e — > 1/e+const. The shape function, in contrast 
to naive expectations, is not a physical distribution, but it is affected by 
regularization scheme effects even at the leading, double-logarithmic level. 
We show, however, that it factorizes most of the non perturbative effects in 
a class of regularization schemes. 

This paper is devoted to a wide audience, i.e. not only to perturbative 
QCD experts, but also to phenomenologists who are interested in the field 
theoretic aspects of this area of B physics, as well as to lattice-QCD physicists 
who may wonder about the possibility of simulating the shape function. We 
have therefore tried to give a plain presentation of our method, together with 
a self-consistent description of the known results to be found in the literature. 

In section [2] we give a simple introduction to the physics of semi-inclusive 
heavy flavour decays. In section |3] we present our strategy, based on factor- 
ization, in order to consistently combine perturbative and non-perturbative 
contributions and to arrive at a formal definition of the shape function in 
field theory; we outline the main steps and the relevant issues. In section § 
we review the standard OPE derivation of the shape function in the effective 
theory; this section can be skipped by the experienced reader. In section |5] 
we return to the strategy outlined in section |] and apply our factorization 
procedure in the quantum theory to a specific class of loop corrections. Our 
technique is completely general, but we believe that it is better illustrated 
by treating in detail a simple computation, which illustrates most of the 
general features. In section |6] we discuss factorization in the framework of 
the effective theory on the light-cone, the so called Large Energy Effective 
Theory (LEET), which is the relevant effective theory for these processes at 
low energy. In section [7] we describe the properties of the shape function in 
the effective theory in various regularizations and its evolution with the UV 
cutoff or renormalization scale. We also discuss our results on factorization 
and clarify a controversial factor of 2 in the evolution kernel of the shape 
function. Section |] contains the conclusions. 
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2 Physics of semi-inclusive heavy flavour de- 
cays 

Let us begin by discussing Fermi motion. This phenomenon, originally dis- 
covered in nuclear physics, is classically described as a small oscillatory mo- 
tion of the heavy quark inside the hadron, due to the interaction with the 
valence quark; in the quantum theory it is also the virtuality of the heavy 
flavour that matters. Generally, as the mass of the heavy flavour becomes 
large, i.e. as we take the limit we expect that the heavy particle decou- 
ples from the light degrees of freedom and becomes "frozen" with respect to 
strong interactions. That is indeed true in the "bulk" of the phase space of 
the decay products, but it is untrue close to kinematical boundaries, as in re- 
gion ([I]). This is because a kinematical amplification effect occurs, according 
to which a small virtuality of the heavy flavour in the initial state produces 
relatively large variations of the fragmentation mass in the final state. To see 
how this works in detail, let us begin with a picture of the initial bound state. 
We assume that the momentum exchanges r M between the heavy flavour and 
the light degrees of freedom are of the order of the hadronic scale, 

|r„| ^0(A QC d), (12) 

as we take the infinite mass limit ([5]). In other words, we assume that the 
momentum transfer does not scale with the heavy mass but remains essen- 
tially constant^. This assumption, which is rather reasonable from a physical 
viewpoint, is at the basis of the application of the HQET fl4}l . Let us discuss 



for example the decay (0). The initial meson has momentum 

Pb = m B v, (13) 

where v is the 4-velocity, which we can take at rest without any loss of 
generality: v ^ = (1; 0, 0, 0). The final hadronic state X has a momentum 



and invariant mass 



Q = m B v - q (14) 



m\ = Q 2 . (15) 



3 We must specify that we consider an initial hadron containing a single heavy quark: 
hadrons containing more than one heavy quark, such as for example quarkonium states, 
need a different theoretical treatment [O) . 
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In eq. flH| ) g M is the momentum of the virtual W or, equivalently, of the 
leptonic pair. We isolate in the decay a hard subprocess consisting in the 
fragmentation of the heavy quark. If the valence quark - in general the light 
degrees of freedom in the hadron - have momentum —k', the heavy quark 
has a momentum^ 

Pq = m B v + k' (16) 

and a virtuality 



m 



B 



2m>BV ■ k' + k' 7^ (in general) . (17) 



The final invariant mass of the hard subprocess, i.e. the fragmentation mass, 
is 

m\ = (p Q - qf = (Q + k'f = m 2 x + 2Q-k' + k' 2 ~m 2 x + 2Q-k' (18) 

and this is the mass that controls the kinematic of the hard subprocess, i.e. 
the Sudakov form factor (the difference between m 2 x and m 2 x is that we do 
not include in the latter the momentum of the valence quark). The term 
k' 2 has been neglected in the last member of eq. (18) because it is small, 
as gluon exchanges are soft according to the assumption flT^). We take the 
motion of the final up quark in the —z direction, so that the vector Q has 
large zero and third components, both of order E x , and a small square; we 
have therefore for the average in the meson state: 

(Q-k') = Q-(k')~0(E x A QCD ). (19) 

A fluctuation in the heavy quark momentum of order Aqcd i n the initial 
state produces a variation of the final invariant mass of the hard subprocess 
of order 

5m 2 x ~ O (Aqcd E x ) ■ (20) 

An amplification by a factor E x has occurred, as anticipated. The fluctuation 
Q20D is of the order of (0) and so it must be taken into account. 

We will discuss the shape function at length in sections f| and |7|, but let 
us introduce now some of its more important properties. If we consider a 



For the appearance of tub instead of mb, see footnote in section 4.1 
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heavy quark with the given off-shell momentum ([16]), we find for the shape 
functionF] 

f(k + ) part = 5 (k + - k' + ) + O (a s ) , (21) 

where 

* + ^ (22) 

Selecting the hadronic final state, i.e. k+, we select the light-cone virtuality 
k' + = k + of the heavy flavour which participates in the decay. After inclusion 
of the radiative corrections, we find that in general k' + > k + Equation (PT|) 
is analogous to the relation between the Bjorken variable xb = —q 2 / (2p • q) 
(p is the momentum of the hadron and q that of the space-like photon) and 
the momentum fraction x of partons in the naive parton model, where we 
have 

q(x) part = 5(x-x B ) + 0(a s ). (23) 

In this case, as is well known, by selecting final state kinematics, i.e. ig, 
one selects the momentum fraction x = xb of the partons that participate in 
the hard scattering. Just as in the heavy flavour decay, radiative corrections 
lead to a softening of the above condition in x > xb, due to the emission of 
collinear partons. 

We note that even with the amplification effect fl2~0|), Fermi motion effects 
are irrelevant in most of the phase space, where typical values for the final 
hadron mass are 

m 2 x ~0(E x ). (24) 

This is in agreement with physical intuition. 

As will be proved in section |4.4j , the shape function can be interpreted as 
the distribution of a variable mass. The virtuality of the heavy flavour can 
be represented by a shift of its mass, mb — > m*. In other words, an off-shell 
particle with a given mass, i.e. with the momentum flTBD , can be replaced by 
an on-shell particle with a variable, virtuality dependent, mass, i.e. with a 
momentum m* (k + ) v. The physical distribution is obtained by convoluting 

5 The final state consists of a massless on-shell quark at the tree level. 

6 To obtain the hadronic shape function, the "elementary" or "partonic" shape function 
in eq. ( ^l| ) has to be convoluted with the distribution ip (fc+) of the primordial light-cone 
virtuality k' + of the heavy quark inside the hadron. 
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the distribution of an isolated quark of mass m* with the probability distri- 



bution for such a mass (see eq. fl74D): this is the basis of the factorization 
theory for the semi-inclusive heavy flavour decays. 

Fermi motion is a non-perturbative effect in QCD because it involves low 
momentum transfers to the heavy flavour (cf. eq. (|12"D), at which the coupling 
is large; it does however occur also in QED bound states, where it can be 
treated with perturbation theory []. 

The second phenomenon relevant in region ([]]) is related to soft gluon 
emission and it is of a perturbative nature - it is a case of the Sudakov form 
factor in QCD fI5|j . The up quark emitted by the fragmentation of the heavy 
flavour with a large virtuality - of the order of the final hadronic energy Ex 
- evolves in the final state, emitting soft and collinear partons, either real 
or virtual. Since the final state is selected to have a small invariant mass 
(cf. eq. (^)), real radiation is inhibited with respect to the virtual one. That 
means that infrared (IR) singularities coming from real and virtual diagrams 
still cancel, but leave a large residual effect in the form of large logarithms^: 



log m x I E 



m 2 x /E 



<** 17 .2 • (26) 



Schematically, the rate for final states with an invariant mass m x has double- 
logarithmic contributions at order as, of the form: 

f Ex f 1 ded9 2 c / 2 m 2 \ 
real = a s / / -— 5 e9 2 - (27) 



o Jo e # 2 V E x 



and 



Jm 2 \ f Ex [ l dedd 2 

where e is the gluon energy, 9 is the angle between the up and the gluon, and 
= 7r — 9 is the polar angle of the gluon 3-momentum. The perturbative 



7 Consider for instance an atom composed of a /i and an e, decaying by /i fragmentation. 
8 The plus-distribution is defined as 

l0g ^ =6{ X ) l ^-5{ X )[ ll ^dy. (25) 

+ x Jo y 
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corrections of the form ( p6|) blow up at the Born kinematics mx = 0, which is 
the threshold of the inelastic channels. For this reason, the above corrections 
are often called threshold logarithms and need a resummation to any order 
in as- 

3 Overview of Factorization 

The aim of this paper is a detailed study of factorization in semi-inclusive 
heavy flavour decays and of the properties of the shape function in field the- 
ory In order to trace all the perturbative and non-perturbative contributions 
to the process, it is convenient to perform the factorization in two steps. In 
the first step the heavy flavour is replaced by a static quark. That is ac- 
complished by taking the infinite mass limit (|5]), keeping Ex and mx fixed. 
With this, the hadronic tensor loses a kinematical scale, namely the heavy 
flavour mass uiq: 

W^ v {m Q , E x , m x ) -> W^ u (E x , m x ) , (29) 
where the effective hadronic tensor is defined as 

W, v = Y,(H Q \Jl\X)(X\J,\H Q )5\p B -q- P x), (30) 
x 

and it contains the static-to-light currents 

J^x) = q(x)T^Q(x). (31) 

The difference between the two tensors in eq. (^) is incorporated into a 
first coefficient function or hard factor. While in full QCD the vector and 
axial currents are conserved, or partially conserved, so the renormalization 
constants are UV-finite and anomalous dimensions vanish, this property does 
not hold anymore in the HQET: the effective current with a static quark is 
not conserved and it acquires an anomalous dimension 7} f][[L6]: 



(/ +f J )(^) = 0. (32) 



dlogjj 



9 In eq. (^2|) and (|3^), we are representing the evolution schematically, without details; 
f.i., we do not distinguish between the anomalous dimensions of the vector and axial 
currents. 
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As a consequence also the hadronic tensor acquires an anomalous dimension, 
which equals twice that of the vector or axial current: 



dlog/i 



+ 2fj = 0. (33) 



All this is very easily understood by observing that the original QCD tensor 
W^y is UV-finite at one loop but it does contain as log rriQ terms, and so 
it is divergent in the infinite- mass limit (|5]). If this limit is taken ab initio, 
i.e. before regularization, these terms manifest themselves as new ultraviolet 
divergences, an heritage of the log mq terms of the original tensor. We may 
say that the dependence on the heavy mass is promoted to UV divergence; 
in practice 

a s log — >a s \og— , (34) 

where Ai is an UV cutoff if we deal with the bare theory, or a renormalization 
point if we deal with the renormalized theory; in principle Ai <C mq. At the 
end of the game, the effective hadronic tensor still depends on three scales, 
just like the original one, 

W, u = W flu (E x ,m x ;A 1 ). (35) 
The original tensor Wu U is parametrized in terms of five independent form 



factors ||17|| . For the HQET hadronic tensor (|30|) there are instead relations 
between the form factors originating from the spin-symmetry of the HQET. 
In particular, the structure in logmg/ Ex of the original QCD tensor can be 
understood by looking at the UV divergences of W^ u @. 

After the first step W^ v still contains perturbative contributions. The 
latter are factorized with a second step, which corresponds to the limit (^). 
Additional UV divergences are introduced also with this second step, which 
must be regulated with a new cutoff A 2 . In principle A 2 <C Ex, since Ex — > 
oo. As before with the heavy mass logarithms, soft and collinear logarithms 
are promoted to ultraviolet logarithms: 

/logro|/B|\ /logHW^ 

n^7irJ + H -w* ) + - ' ' 



J In Dimensional Regularization (DR), this means simple poles 1/e. 
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The second factorization step involves double-logarithmic effects of an in- 
frared nature, in contrast with the single logarithms of the large mass of 
the first step. In practice, we separate perturbative contributions from non- 
perturbative ones starting with a cutoff 

A 2 ~ E x (37) 

and lowering it to a much smaller valueQ 

A' 2 < E x . (38) 

The contributions of the fluctuations with energy between A 2 and A' 2 are 
put into a second coefficient function, while the contributions below A' 2 are 
factorized inside the shape function. The latter is defined in the framework 
of a low-energy effective theory, with a cutoff given by 

A ET = A' 2 . (39) 

Most of the non-perturbative effects in lattice-like regularizations are con- 
tained in the shape function, which uniquely determines the final, non- 
perturbative, hadronic tensor 

W, u ^^(F Q |JjX)(X|J M |if Q )5 4 (p B -g-p x ), (40) 

x 

containing the effective-heavy-to-effective-light currents 

J f ,(x) = q(x)T l _ l Q(x). (41) 

It is worth noting that the tensor ( [4"0"D involves a single form factor, propor- 
tional to the shape function itself (see eq. (|65|)). That is again a consequence 
of the spin-symmetry of both HQET and LEET jrj|, which is more effi- 
cient than that one of the HQET alone. The shape function is completely 
non-perturbative and perturbative factors can no longer be extracted. 

The effect of lowering the UV cutoff (eqs. ( p?|) and (^)) is incorporated 
inside a coefficient function, which, unlike more simple cases such as the light- 
cone expansion in DIS, is not completely short-distance dominated. Some 

n In order to avoid substantial finite cutoff effects, the condition A' 2 3> Aqcd must hold. 
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long-distance effects are left in the coefficient function, but they are expected 
to be suppressed on physical grounds. Finally, the introduction of ultravio- 
let divergences with factorization, implies scheme- dependence issues for the 
shape function, which are rather dramatic because of the double-logarithmic 
nature of the problem (cf. eqs. (p7|) and (|£f)). 

In fig. 1, we give a pictorial description of the above procedure. 



4 OPE 

The amplitude for the decay (|7|), which we take as our example from now 
on, can be written at the lowest order in the weak coupling as 

A = ^(lv\L li \0)(X\J»\B), (42) 

where L M is the leptonic current and J M is the hadronic one: 

J^x) = q(x)T^b(x) (43) 

with r M = 7 M (1 — 75), q(x) a light quark field and b(x) the beauty quark 
field. Taking the square of (|4*2|) and summing over the final states, we arrive 
at the hadronic tensor defined in eq. (§). By the optical theorem, we can 
relate the hadronic tensor to the imaginary partQ of the Green function 
or forward hadronic tensor T^ v : 

= --Im TV, (44) 

7T 

T, v =-ij d 4 xe-^(B\T (Jl(x)J„{0)) \B). (45) 



where 



4.1 The HQET 

We are interested in the evaluation of in the effective theory and we 
discuss in this section the first factorization step: replacing the beauty quark 



12 Since T M is in general complex, we should say, more properly, the absorptive part. 
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by a quark in the HQET. As is well known, we can decompose the heavy 
quark field b(x) into two effective quark and antiquark fields h v and H v 

b(x) = e - imBV - x [h v (x) + H v (x)], (46) 

satisfying 

P + h v = h v , P. h v = 0, P_ H v = H v , P + H v = 0, (47) 

where P± = (1 ± v)/2 are the projectors over the components with positive 
and negative energies, respectively. The field H v is neglected (which amounts 
to neglecting heavy-pair creation), so that 

b(x) ~ e - imBV - x h v (x). (48) 

By using eq. ( pEB] ) we obtain 

%v = ~i\ d*x e iQ - x (B(v)\T h v (x)Tl q(x) g(0)r, h v (0)\B(v)). (49) 



We now use the Wick theorem and we single out the only contraction that 
is relevant to B decay: 

Tp, = J d A x e iQx (B\h v (x)TlS(x\0)T u h v (0)\B), (50) 

where S(x\0) is the light quark propagator. Note that the operator entering 
the right-hand side of eq. (|50D is already normal ordered, since h v has only 
the component that annihilates heavy quarks, while h v only the components 
that create them. We can express the Fourier transform of the light quark 
propagator as0 

S(Q + iD) = = _^ /2 ___ _ _ (51) 



13 We prefer to refer to the physical i?-meson mass rather than to the unphysical 6-quark 
mass. Their difference is of order Aqcd, so it is a 1/tob correction and can be neglected 
in our leading-order analysis. Furthermore, in perturbation theory there is no binding 
energy so that mj, = nis- 

14 The notation is very compact. For more explicit representations of the propagator see 
ref. H. 
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where = i/2[ r y fJi , 7„] is a generator of the Lorentz group, G^ = —i/g[D^, D v ] 
is the field strength and = <9 M — is the covariant derivative. In eq. 
(|5"TD, denotes, as usual, an infinitesimal positive number and gives the 



prescription to deal with pole or branch-cut singularities. There are three 
different regions according to the value of the jet invariant mass, which are 
described by three different full or effective theories: 

i) m\ ~ Q{E\\ 



ii) m\ ~ 0{A QCDh 



in) m 



x 



2 

Q< 

0(A QCD E x ). (52) 



Since the derivative of the rescaled h v field brings down the residual mo- 
mentum k', and it is therefore an operator with matrix elements of order 
O(Aqcd), the matrix elements of the operators entering the light quark 
propagator have a size of the order of 

{iD) ~ O(Aqcd), 

(2W-Q) ~ 0(A QCD E x ), 

(D 2 ) ~ 0(Al CD ), 

K„G"-> ~ 0(AJ CD ). (53) 

Let us discuss these regions in turn in the next section. 



4.2 General kinematical regions 

i) This region corresponds to a jet X with a large invariant mass, of the 
order of the energy: 

m x ~ O (E x ) . (54) 

To a first approximation all the covariant derivative terms can be ne- 
glected, so that 

S{Q + ' D) ^ Q^+iO' (55) 

i.e. the light quark can be described as a free quark. A higher accuracy 
is reached when expanding the propagator in powers of the covariant 



14 



derivative operators up to the required order. We have here an applica- 
tion of the l/m,B expansion up to a prescribed (finite) order [^]. In this 
region there are no large adimensional ratios of scales, the latter being 
all of the same order. This implies that in perturbation theory we do 
not hit large logarithmic corrections to be resummed to all orders in as- 
This region is not relevant to the endpoint electron spectrum because 
the hadronic jets takes away most of the available energy. This region 
will not be discussed further here. 

ii) This region involves a recoiling hadronic system with a mass of the order 
of the hadronic one: it can be a single hadron or very few hadrons. 
The dynamics is dominated by the emission, with consequent decay, of 
few resonances; it is a completely non-perturbative problem. Accord- 
ing to the estimates (|53|), no term can be neglected in the light quark 
propagator. We are faced with full QCD dynamics as far as the final 
hadronic state is concerned. This region must be evaluated by an ex- 
plicit sum over all the kinematically possible hadronic states, and the 
latter have to be computed with a non-perturbative technique such as 
a quark model or lattice QCD. This region will not be discussed here 
either. 

Hi) This region is intermediate between i) and ii) and as such it has both 
perturbative and non-perturbative components. Roughly speaking, we 
have to take into account non-perturbative effects for the initial state 
hadron, while we can neglect final state binding effects. This region is 
characterized by a small ratio of the jet invariant mass to the jet energy, 
and thus involves the large adimensional ratio in (|3|). As always is the 
case, perturbation theory generates logarithms of the above adimen- 
sional ratio, eq. (p6|). The term 2iD ■ Q at the denominator cannot be 
brought at the numerator (with a truncated operatorial expansion) be- 
cause it is of the same order as m 2 x . At lowest order, the other covariant 
derivative terms can be neglected, to give: 

S(Q + iD)~— 9 Q-— . (56) 

vv ' m 2 x + 2iD-Q + iO K ' 



15 It is clear that a consistent inclusion of the 1/ms corrections involves also the expan- 
sion of the heavy quark field b(x) into the effective quark field h v (x) up to the required 
order. 
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One can reach a higher level of accurary keeping these latter corrections 
up to a given ordeij^]. The rest of the paper deals with region Hi) at 
the lowest order in l/m B . 



4.3 The LEET 

In this section we discuss the second factorization step, which involves the 
description of the final up quark in the LEET, according to eq. fl5"6"|). Let us 
define the adimensional vector 

< 67 > 

This n M has a normalized time component, n = 1. In the "semi-inclusive" 
endpoint region Hi): 

= o(^jgr) «i- (58) 

We will show later that n can be replaced by a vector lying exactly on the 
light-cone, i.e. 

n — > n, (59) 

where = (1; 0, 0, —1) (n 2 = 0), representing the direction of the hadronic 
jet, the —z axis. We can write 



where k + has been defined in eq. (|22| ) and D + = n ■ D. We can simplify the 
tensor structure of T^ u by using the identity 

h v T^h v = ^Tr(r M P + ) h v h v - ^Tt(j^ 5 P + T ^P + ) h^^hy, (61) 



16 We envisage a relation between the 1/Ex corrections to the shape function and the 
power-suppressed perturbative corrections of the form as /Ex log 2 (Ex/mx)- 
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which is valid for any T^. The matrix element of the axial vector current 
between the £>-meson states vanishes by parity invariance, so that 

?V = v F(k+), (62) 

where we have defined the "light-cone" function 

F(k + ) = (B{v) | ht iD+ _\ + + i0 hv I B(v)), (63) 



TlQT u P + (64) 



and 

s ^ = 2 Tr 

is the "spin factor", containg the leading spin effects. The factor l/(2v ■ Q) 
is a jacobian, which appears as we go from the full QCD variable Q 2 to the 
effective theory variable k + . 

Taking the imaginary part of T^, we obtain (see relation (f£|)): 

= s ^ f( k +)i (65) 

where 

f(k + ) = --lmF(k+) (66) 

7T 

is the shape function. By using the formula 
1 „ 1 



iD + — k+ + iO i D + — k. 



iir5(iD + - k+), (67) 



we recover the definition of the shape function given by eq. (|10]). Note that 
it involves the non-local operator hi S(k + — iD + ) h v , which results from the 
resummation of the towers of operators of the form (Q ■ D) n . 

17 A physical argument for the spin factorization is that, in the limit mg — > oo, the 
spin interaction of the 6-quark in the B meson vanishes; therefore we can average over the 
hclicity states of the b quark [M . 
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4.4 The Variable Mass 



The hadronic tensor can be written in the effective theory in terms of the 
shape function as: 



dk + 5 (Q 2 + 2k + v ■ Q) /(*+); 



(68) 



in the second member, k + is an integration, i.e. dummy, variable. In the free 
theory, with an on-shell 6-quark (i.e. k! = in eq. (|16"D), 



so that 



f°{k+)=S{k + -0), 



K> = v W 2 -o). 



(69) 



(70) 



The hadronic tensor can be written, up to terms of order k\ ~ O (Aq C£) ) 



as 



W 



[IV 



Sfiu(Q) I dm* 5 (Ql - 0) /(m* - m B ) 



(71) 
(72) 



where we have defined 

Q* = m*i> - q 

and m* is the "variable" or "fragmentation" mass, defined in eq.(ll]). Since 
m* is just a shift of k + , the range is 



(73) 



Inside we can replace Q with Q*, because that amounts only to a correc- 
tion of order k + = 0(Aq C d), so thatQ 



m B 



dm* (p(m*) W°(v, Q* 



(74) 



where 



<>,Q*)= s^(Q*)6(Ql-0) (75) 
is the hadronic tensor in the free theory for a heavy quark of mass m* and 



(p(m*) = /(m* - m B ) 



(76) 



18 We replace by the lower limit of integration, because the relevant region is m* 
m B - O(Aqcd)- 
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is the distribution for the effective mass m* of the 6-quark inside the B- 
meson at disintegration time. Equation (74) is the fundamental result of 
factorization in semi-inclusive heavy flavour decays: it says that the hadronic 
tensor in the effective theory can be expressed as the convolution of the 
hadronic tensor in the free theory with a variable mass times a distribution 
probability for this mass. That offers also the physical interpretation to the 
shape function anticipated in the introduction: it represents the probability 
that the b quark has an effective mass m* at the decay time. Since this 
tensor encodes all the hadron dynamics, any distribution can be expressed 
in a similar factorized form. 



5 Factorization in the quantum theory 

In this section we discuss factorization in the quantum theory, i.e. the separa- 
tion of short-distance and long-distance contributions, including loop effects. 

A shape function f(k + )® CD and a light-cone function F(k + )® CD can also 
be defined in full QCD by means of the relations [[l(J : 



and 



where As^ and As' are defined as the part of the spin structure not pro- 
portional to s^ u . The tensors As M „ and As^, represent residual spin effects 
not described by the effective theory (ET), which do not contribute to the 
Double-Logarithmic Approximation (DLA)Q. In DLA the forward tensor can 
therefore be written as 

T* CD = V ^ F(k + f CD , (77) 

where the "light-cone function" is given by 

F(k + f CD ee -^1— [1 + aC] , (78) 



in 



Note that As M „ and As' are, in general, different; this was not noted in |lC 
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a = asCp/ir and C is the scalar triangle diagram (see fig. 2): 



_• O / — 1 1 1_ 

We have set the light quark mass equal to zero |]2T| . 

The hadronic tensor relevant to the decay is obtained by taking the imagi- 
nary part according to eq. (fHj). This transforms the products in convolutions, 
which are converted again into ordinary products by the well-known Mellin 
transform 122] . 



Infrared singularities (soft & collinear) are regulated by the virtuality 
Q 2 of the external up quarkF]. We may write 



2 2 \ / 2 

n „ „ n \ „ / n 



Q^ = E X (1 + T ;0 ) 0,-1 + T j = £ x + (80) 



and 



n 2 



n» = v» + —v% , n» = vt (81) 



where we defined the light-cone versors 

v + = (l; 0,0,1), u_ = (l; 0,0,-1). (82) 

Let us now consider the properties of the integral C. First, it is adimen- 
sional. Second, it is UV-finite for power counting: the integrand has three 
ordinary scalar propagators with a total of six powers of momentum at the 
denominator. This implies that C does not depend on an ultraviolet cutoff 
Auv as long as it is larger than any physical scale of the process, namely 

Auv > m B . (83) 

Third, as already discussed, C is also IR-finite as long as Q 2 ^ 0. For Q 2 > 
there is an imaginary part, related to the propagation of the real up and 
gluon pair, while for Q 2 < the integral is real. Therefore C does depend 
on adimensional ratios of three different scales: rriB,Ex and mj. There are 
only two independent ratios, which we choose as ms/Ex and nix /Ex- We 

20 This is consistent because a virtual massless quark is not degenerate with a quark and 
a soft and/or collinear gluon. 
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are going to decompose the integral C in a sum of various integrals; at the 
end, one of them will correspond to the double-logarithmic contribution to 
the shape function f(k + ) in the low-energy ET. The other integrals represent 
additional contributions and they are mostly short-distance dominated in 
lattice-like regularizations. This decomposition consists of two separate steps, 
which will be described in the following sections. 



5.1 From QCD to HQET 

In the first step we isolate a hard factor by simply subtracting and adding 
back the integral with the full beauty quark propagator replaced by a static 
one (see fig. 3) 

C = C s + C h , (84) 

where 4 

C s = -ivQj— {l + Q )2+iQ ^yV^ PTM> (85) 

and 

C ~i V ' Q [ dH 1 - 1 (86) 

h % 2m J vr 2 (l + Q) 2 + iO vl + iO v ■ I + I 2 /2m + iO' 1 ' 

The above decomposition parallels that one performed in section [| in oper- 
atorial language. The light-cone function factorizes at order as according 
to 

We expect that C s has the same infrared behaviour as C; it will be subjected 
to a further decomposition in the next sections; Ch is the "hard factor", 
i.e. the difference between QCD and the static limit for the b quark. The 
latter integral is both UV- and IR-finite. The ultraviolet finiteness stems 
from power counting: the integrand has two scalar propagators and a static 
propagator with a total of five powers of momentum in the denominator PJ 
In the infrared region, all the components of the loop momentum are small 

IRj_ l»^pl» P^O, (88) 

21 It is known that ultraviolet power counting may fail in effective theory integrals 
when there are LEET propagators because of the occurrence of an ultraviolet collinear 
region p0|. The integral Ch, however, contains only an HQET propagator. 
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so we can neglect the terms that are quadratic in Z M in the propagator de- 
nominators: 



Ch ' IR ~ J dl 2l-Q + Q* + iO (vl + ioy {89) 

Integrating over Iq, and closing the contour in the upper half of the /o-plane, 
we see that there are no enclosed poles and the integral vanishes (QED). 
Inside Ch we can therefore make the replacement Pi 

Qf _> Q^ = E x v^ (g 2 = 0). (90) 

It follows that Ch depends only on m B and E X - Ch = Ch{rriB, E x )- Since it 
is adimensional, it may depend only on the adimensional hadronic energy 



2E, 



(91) 



m B 

i.e. Ch = Cji(z). An explicit computation gives 

C h = log(l - z) log z + Li 2 {z) ~ -zlogz (z<l), (92) 

C/i does not contain large logarithmic contributions in the limit z — > 0, i.e. 
log z terms 0. This is related to the fact that Ch and C s are UV- convergent. 
In general, single logarithms of the hadronic energy, log 2, do appear, rep- 
resenting the difference between the interaction of a full propagating heavy 
quark, of mass m#, and that one of a static quark. The relevant interaction 
energies are between the beauty mass m B and the hadronic energy Ex, 

[ m B dk 2 

as — = -2a s log z. (94) 



22 With this substitution, terms related to higher twist contributions of the forms 
(mx /ttlb)" 1 and (mx / Ex)™ 2 are neglected, but this is in agreement with our leading-twist 
ideology (the indices n\ and rii are integers). 

23 The dilogarithm is defined as 

Ln{ z )^-f l ^pl^ = ±^ (|*| <1). (93) 
Jo ^ n 
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The logarithms are resummed, as usual, by replacing the bare coupling 
with the running coupling and exponentiating, so that the above formula is 
corrected into: 



= exp [-2 ct 5 7o log z + 27 /3 c? s log 2 z H ] ,(95) 

where j3 = l/(47r)(ll/3 N c — 2/3 nf) and as = as(mfl). 

Let us summarize the above discussion. A first coefficient function is 
introduced, which takes into account the fluctuations with energy e in the 
range 



In the language of Wilson's renormalization group, we are lowering the UV 
cutoff of the effective hamiltonian from to Ex- 

5.2 Infrared factorization 

The second factorization step involves the separation of the various infrared 
contributions to the process, one of which will ultimately lead to the shape 
function. This step forces us to introduce explicitly an ultraviolet cutoff 
A from which the various factors depend separately. In other words, the 
decomposition of C s introduces fictitious UV divergences, which cancel in 
the sum. As anticipated in the overview section, infrared factorization in our 
scheme involves two different operations: 

• the separation of the various pole contributions to the QCD amplitude 
according to the Cauchy theorem; 

• the lowering of the UV cutoff from Aw ^ Ex to Auv = ^et <C 
Ex, where Aet is the cutoff of the final low-energy effective theory, in 
which the shape function is defined (the latter contains the majority of 
the non-perturbative effects). 

Step 1) will be discussed in this section, while step 2) is treated mostly 
in the next section. 




E x < £ < rn B . 



(96) 
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C s is UV-convergent, as is clear again from power counting, and it does 
not depend on the beauty mass tub, which has been sent to infinity, so that 
C s = C s (Ex,rrix)- The only adimensional variable that can be constructed 
out of Ex and m x is their ratio or, equivalently, n 2 (see eq. (0)). Since 
C s is adimensional it may depend only on n 2 : C s = C s (n 2 ). The explicit 
computation in DLA gives 

C s = ~ \og 2 (-n 2 - iO) (DLA) . (97) 

The infrared factorization is performed by integrating C s over the energy In 
using the Cauchy theorem. There are three poles in the lower half of the 
Z -plane related to the propagation of a real static quark, a real gluon and a 
real up quark, located respectively at 



l = -iO, l = +\l\-iO, l = -Q + ^(Q 3 + l 3 y + ll-i0. (98) 

In the upper half-plane, instead, there are only two poles, related to the gluon 
and the up propagator: 



In = -\l 



+ iQ, l = -Qo-\/(Q3 + h) 2 + ll + tO (99) 



The poles in fl9"5|) are conventionally related to a propagation that goes back- 
ward in time; they will therefore be called the antiparticle poles. We close 
for simplicity the integration contour in the upper half-plane and we have 
two residue contributions related to the antigluon pole ad the anti-wp pole 
respectively (see fig. 4): 

C 8 = C g + C q . (100) 
The antigluon and the anti-wp contributions are given respectively by 



Cg = *v-Q [ d 3 l ' 1 1 



vr ^ J i Q - | i | + v ■ I + iO Q 2 + 21 ■ Q + iO 

2 



i =-| I |+i0 



Cg = -V-Q I d 6 l 



h + Qo-\l + Q\+i0v-l + t0l 2 + i0 



(101) 



lo=-Qo— s/(Q3+h) 2 +l 2 x +iO 



As we see by power counting, C g and C q are UV-divergent and it is therefore 
necessary to introduce an ultraviolet regularization to treat them separately. 



24 



The light-cone function factorizes after this second step as 

F(k + f CD = -1— [1 + a C h ] [1 + a C q ] [1 + a C g ] , (102) 
— k + + 20 

i.e. as a product of three factors. 
5.2.1 Wilson line representation 

Before explicitly computing these 3-dimensional integrals, let us represent 
them as 4-dimensional ones, i.e. as one-loop integrals of a properly chosen 
field theory: 

/cfc^l 1 1 1 

-i Q2 + 2i-Q + io V^TTm p + iO (103) 

and 

/d/^l I 1 I 

^2 (/ + Q)2 + , vl + zO Q 2 + 2l-Q + i0 (1 ° 4) 

The proof of the above equations is just by integration over l : closing the 
integration contour in the upper half-plane of lo, we enclose a single pole, 
whose residue gives the 3-dimensional integrals in eqs. ( |101|) ; C g and C q 
involve one full - i.e. quadratic - propagator and two eikonal - i.e. linear - 
propagators. It is easy to check that the algebraic sum of C g and C q in the 
above expressions reproduces the integral C s defined in eq. (§§). Introducing 
the variable k + , the integral C g can be written in the "familiar" form 

r (h \- ~ l I dAl 1 1 1 

C 9 (k + ) = — J — _ k+ + lmn + iQ v . l + lQ p^Q- (105) 

The geometrical interpretation is the following: C g (k + ) is the one-loop cor- 
rection to a vertex composed of an on-shell Wilson line along the time axis, 
and a Wilson line along the direction n off-shell by k + . Note that 

l-n = l+ + — L, (106) 
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so that eq. ( |105| ) represents the vertex correction in Feynman gauge to the 
function 

F(M„vo - (B(v) I ft t__l__ ft „ | B{v)) . (107 ) 

The imaginary part of the above function equals —1/tt times the shape func- 
tion off the light-cone, n 2 ^ 0: 

f(k + ) n ^ = (B(v) | hi 5 (k + - iD + - in 2 /AD J) h v | B(v)). (108) 

In the limit n 2 — > (see section ||) we recover the correction to the light-cone 
function F(k + ), which was computed in ref. ||, |], [ToU . 

We can give a similar description for C q . The 4-dimensional representa- 
tion for C q involves an on-shell Wilson line along the time direction, a Wilson 
line along the direction n off-shell by k+, and a light quark propagator with 
momentum I + Q 

i [ dH 1 1 1 

q= 2j^ -k+ + l-n + i0 vl + iO (l + Q) 2 + iO' (1 ° 9) 

Note that the expressions for C g and C q are very similar: they differ by an 
overall sign and by the replacement in the quadratic propagator of I — > l + Q. 
For future reference, let us note that the latter shift involves only the zero 
and the third components of /, not the transverse ones. 

In fig. 5 the decomposition of C s into C g and C q is represented. 

5.2.2 Space Momenta Cutoff 



We consider the bare theory with the regularization introduced in ref. [ 10 1 : a 
sharp cutoff on the spatial loop momenta As (see next section). Integrating 
C g over Iq by closing the integration contour upward and integrating over the 
azimuthal angle, we obtain 

Ms rl i 

C„ = - I dl dcos9- ; —— -, (110) 

9 Jo J_x k + + l(l-cos9)+n 2 l/4(l + cos9) , K ' 



where we have used the definition of Q M in eq. (|80|) and I = | I \ . Integrating 
over the polar angle, we obtain 



2A S 



dl 



C 9 = ~l ylog 



o 



k+ + l 



_k + + n 2 l/A_ 



(111) 



26 



Assuming a cutoff much larger than any physical scale in the process, i.e.0 



As^E x , (112) 

we obtain in DLA Q 

r< 1 i 2 E x , A s E x , x 

2 fc + - 2O £x — fc+ + iQ 

Three scales enter in C g \ k + , A s and E X - The appearance of k + and the 
cutoff A s was expected, because these two quantities represent the infrared 
and the ultraviolet scale, respectively. The noticeable fact is that also the 
hadronic energy E x makes its appearance. C g contains a double-logarithm 
of the infrared kinematical scale k + (related to the overlap of the soft and the 
collinear region, which extends up to As); it also contains a single logarithm 
of the cutoff. The appearance of the hadronic energy Ex comes from the 
necessity of a third mass scale for the function C g , which behaves like log 2 k + 
in k + and like log As in As- 

When I 3> Ex the integrand behaves as 



k+ + l 



I l0g [k + + n 2 l/4 



2 



1 , n 

ylog T (H4) 



and produces a single- logarithmic ultraviolet divergence. As eq. (|114j) clearly 
indicates, n 2 ^ regulates the collinear or light-cone singularity: up to now 
we have indeed taken kinematics into account exactly together with a large 
cutoff. It is interesting to note that if we take a cutoff much smaller than the 
hadronic energy (as we will do in the "final" low-energy effective theory), 

Aet < E x , (115) 

we have 

2 

n 2 l < -f Aet < k + , (116) 

24 This is done consistently with the relation (|8^), in which we have taken a large cutoff 
for the computation of C s 

25 

opposite one in C q (see eq. ( |119| ) 



The last member in eq. (|115| ) is an artificial absorptive part that cancels against an 
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and C g simplifies in 

C g (A ET ) ~-J^ - log ( -±— J . (117) 

The quantity n 2 does not enter anymore and the integrand is the same as 
that with the approximate light-cone kinematics n 2 = 0, i.e. with n replaced 
by n. The physical explanation is that soft gluons are not able to distinguish 
between the two slightly different directions n and n. Formally, with the 
small cutoff ( 115| ) we can take the limit (^J) inside the integral. In other 
words, in the low-energy effective theory, we effectively are always in the 
light-cone limit. For I ^> k + , the integrand in eq. ( |117| ) has the asymptotic 
behaviour 

ylog^, (H8) 

implying a double-logarithmic behaviour with respect to Aet upon integra- 
tion over Z, in contrast with the single logarithmic behaviour of the integrand 
in ( 114J) . These properties will be studied systematically in the next section, 
in which we consider the effective theory on the light-cone. 

For the computation of C q , it is convenient to first make the shift I — > l—Q 
in the expression of C q , eq. (|104[) , and then to compute the residue of the 
light quark pole at Iq = — | / | -HO: this is legitimate if condition (|112|) holds. 
We find 

C q = log log — ^— (As^>E x ). (119) 

E x —k+ + ^0 

The three scales appearing in C g do appear also in C q . We note that C q 
contains a single logarithm of k + , i.e. it is subleading by one logarithm in 
the infrared counting with respect to C g . It has a single- logarithmic UV 
divergence. 



If we compute the integral in eq. ( |104| ) with a small cutoff (|115|) , we do 
not find any infrared logarithm, in contrast with what happens instead with 
C g . Thus the logarithmic contributions to C q come from high-energy gluons 
and that is an indication that C q , unlike C g , is short-distance dominated. 

One can check that the correct value of C s is reproduced by summing C g 
and C q ] in particular, UV divergences cancel. 

At the level of logarithmic accuracy, we can replace the strong inequality 
( |112|) with a weaker one: 

As > E x . (120) 
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Setting in particular 

As = E x , (121) 
the expressions for the gluon and quark-pole residue read 

C g (A s = E x ) = -^ 2 j^- = C s , 

C q (A s = E x ) = 0, (122) 

i.e. the gluon-pole term gives the whole contribution while the quark-pole 
factor vanishes. The term C q therefore has the role of correcting C g when 
As ± E x . 

Since ultraviolet singularities are single-logarithmic for a large cutoff (eqs. 
(|113|) and ( |119| )), other regularizations such as DR give similar results. In 



other words, the regularization scheme dependence is the usual one: the 
logarithmic term in the one-loop amplitude is scheme independent while the 
finite part is scheme dependent. 



6 The effective theory on the light-cone, the 
LEET 

In full QCD infrared singularities are regulated by the unique quantity m 2 x ^ 
0. In the effective theory, the light quark propagator is replaced by an eikonal 
one 

1 1 1 



(123) 



(/ + Q) 2 + zO Q2 + 2l-Q + P + i0 Q2 + 2 l-Q + i0' 

In the expression on the right-hand side I 2 has been neglected and as a con- 
sequence enters in two distinct and independent ways: its square Q 2 
represents the virtuality of the eikonal line at = 0, while its components 
Q M are the coefficients of the linear combination of the loop-momentum com- 
ponents 1^ in the term 2 Q^. Unlike full QCD, Q 2 and can be considered 
as independent quantities in the effective theory. We may ask ourselves what 
happens if we take the limit Q 2 — > inside the term 21 ■ Q — > 21 ■ Q while 
keeping Q 2 const ^ 0, i.e. if we make the replacement 

(124) 



Q 2 + 2l-Q + i0 Q2 + 2 l.Q + iO 
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In the usual notation, the above replacement reads 

1 1 



(125) 



— k + + n ■ I + iO —k+ + n ■ I + iO 
corresponding to the limit 

n 2 -»• 0, k + -> const ^ 0. (126) 

This means that we are considering an eikonal propagator that lies exactly 
on the light-cone, with collinear singularities regulated now by k + ^ only, 
instead of by n 2 ^ [0, As we saw before, the limit n 2 — > 

is "invisible" with a small cutoff, simply because the integrand does not 
depend on n 2 in this case. We now want to see what happens in the limit 
( 126 ) with a large cutoff. To perform IR factorization in the light-cone case, 
it is convenient to start from the original QCD amplitude C s in which we 
make the replacement 

(127) 



Q 2 + 2l-Q + P + i0 Q2 + 2/ .Q + /2 + i0 ' 
to obtain 

C s -+C s = -ivQ I ^ 1 r ~pr~ — • (128) 

^ J 7T 2 Q2 + 2 l.Q + p + i o vl + iO P + iO v ; 

It is straightforward to check that C s and C s coincide at the DLA level, i.e. 
that the approximation ( |127| ) preserves the infrared structure, 

C s = C s = -Uog 2 (-n 2 -iO) (DLA). (129) 

The gluon and the quark pole contributions are given in the light-cone limit 
by 



Cg,n 2 =0 = W-Q —T 



dH 1 

tt 2 Q 2 + 2l-Q + i0 v-l + iO P + iO 
dH 1 1 



C ^o - -ivQ / n2 ,^,, -y— H j^-K (130) 



7T 2 Q2 + 2l-Q + P + i0 v-l + iO Qi + 21-Q + iO 
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The above terms are usually called "soft" and "jet" (or "collinear") factor 
respectively, even though we believe that this terminology can be rather 
misleading, as the redistribution of double logarithmic contributions in C g 
and C q is substantially dependent on the regularization. We will show later 
that it is possible, within a specific class of regularization schemes, to confine 
all the double logarithmic contributions in C g . It is immediate to check that 
the two above integrands sum up to the integrand of C s . Making the shift 
I — > I — Q, the quark factor can also be written as 

dH 1 1 1 



C q , n * =0 -ivQj ^ p + Q2 + iQ v , l _ v _Q + iQ Q 2 + 2l .Q + i0 - 

(131) 

6.1 Regularization Effects 

The decomposition of C s into C g n 2 =0 and C q:Tl 2 =0 is strongly dependent on 
the regularization scheme adopted, as a consequence of the fact that double- 
infrared logarithms are promoted to double-ultraviolet logarithms with the 
splitting. We will see that there are substantial regularization scheme effects, 
even for the leading DLA terms. Two different classes of regularizations are 
considered. To the first class belongs the regularization considered in ref. [fL0|1 : 
a sharp cutoff on the spatial loop momenta 

\t\ < As, 

and a loop energy on the entire real axis, 

- oo < l < oo. (132) 

That means, roughly speaking, a discrete space and a continuous time. We 
believe that this regularization gives the same double-logarithm as the ordi- 



nary lattice regularization - the Wilson action [23]. In the latter case all the 



components of the loop 4-momentum are cutoff, not only the spatial ones 

H < A 4 = -, (133) 
a 

where a is the lattice spacing. The physical reason for the equality of the 
double-logarithmic coefficients in the regularizations ( |132| ) and ( |133j ) is the 
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following. Soft and collinear logarithms are both related to quasi-real gluon 
configurations, for which 

k ~ \~T\- (134) 

Cutting off the spatial momenta therefore should cut off also the relevant 
energies as far as soft and collinear singularities are concerned. 

As a representative of the second class of UV regularizations, consider a 
sharp cutoff on the transverse momenta (the x-y plane): 

\ I < Aj_, while — oo < l + , /_ < oo. (135) 

This regularization is "effective" , i.e. it is sufficient to cut on the transverse 
momenta to render the integrals finite. To this class of regularizations belongs 
the Dimensional Regularization (DR), in which most of the effective field 
theory computations have been performed. Let us treat the two cases in 
turn. 



6.1.1 Space Momenta Cutoff 

An explicit computation of the gluon and the quark pole contributions on 
the light-cone in the A^-regularization gives 

r> 1, 2 A s 



2 ° fc + -iO 

C „2-n = -log log . 136 

9 ' n _0 2 B k + -iO 2 B k + -iO v ; 

The behaviour with respect to k+ is the same as in the case n 2 ^ 0. Ultraviolet 
divergences are now more severe than in the case n 2 ^ 0, being of double- 
logarithmic kind. However, the sum is again the correct one 

Cg,n 2 =0 + C 9 ,n 2 =0 = C s . (137) 

In other words, the transition to the light-cone theory implies a rearrange- 
ment of the ultraviolet structure, but the physical observable, C s , is un- 
changed. 
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6.1.2 Transverse momenta cutoff 



The factor C g is better computed in this case by introducing light-cone co- 
ordinates: 

l + = l + l 3 , l_ = l -l 3 . (138) 

Integrating over /_ by closing the integration contour upward and over the 
transverse momentum, we obtain 



C a 



1 



dyy 

1 + y 2 l + nV/4 



log 



1 + 



A 



i + 



n 2 y 2 



Performing the final integration in the case n 2 = 0, we obtain 



a 



g,n 2 =0 



(A ± ) = -log 2 



A 



k + -i0 

For the quark-pole factor C q , the integration over /_ gives 

"A 2 



a, 







rfZI 



1 







x 2 l\ + 2E x x + l k + x+l-i0' 
Integrating over l± we obtain, in the light-cone limit: 



C n 



dx 



1 



withF] 



o x 1 — n 2 x/A — iO 

A, 



log 



A 2 ± x 2 + x + l 
1 + x 



A, 



2£ 



The above integral has two double- logarithmic regions for Aj_ 3> Ex, 

4 1 

(1) : 1< x < — , (2) : =- < x < 1. 

vr Aj_ 

Performing the integration in the two regions, we find 



C g ,n2=0 (A±) = J log 



Ai 



E x (k+ - «0) 



log 2 



E 



x 



26 ^2 



n in the above formula has to be interpreted as — 2k+/Ex- 



(139) 



(140) 



(141) 



(142) 



(143) 



(144) 



(A ± >£ x ). (145) 
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The first double-logarithm on the right-hand side is related to region (1), the 
second one to region (2). For a smaller UV cutoff, we obtain instead: 

C q>n * =Q (A_l) = \ log 2 — ^1— - (E x \k+\ « Al « E x ) . (146) 

Finally, for Aq CD <C <C Ex \k+\, the integral C q vanishes in DLA. 

6.1.3 Comments 



Let us comment on the results ( |14(J| ) and (|145 ) . As with the 3-momentum 
regularization, C g and C q have double-logarithmic UV divergences, again 
a consequence of the light-cone limit n 2 = 0. The most important point, 
however, is that C g has an additional factor of 2 with respect to the spatial 
cutoff case in the coefficient of the double-logarithm of the infrared scale, 
log 2 k + (cf. eqs. ( |113| ) and ( |140| )). With the As regularization, C q has no 



log 2 k+ term, while with the A^ regularization it does. The same double 
logarithm is obtained in the sum C s in both regularizations. In general, the 
appearance of log 2 k + in C q (A±) implies that, with the A^ regularization, 
C q does not describe only collinear contributions but also soft ones We 
interpret this fact by saying that the shape function, in general, does not have 
any physical meaning, but it just represents the gluon-pole contribution to 
a physical process: that result is, as far as we know, new. One generally 
attaches to the shape function a physical meaning - related to the Fermi 
motion; thus, to understand what is happening, we have to start again from 
the beginning. The shape function is obtained from the original QCD tensor 
Wpv considering the infrared limit of small momenta compared with the 
hadronic energy: 

M < E x . (147) 

The tree-level rate in the ET equals the QCD one by construction. However, 
in loops, the condition ( pL47| ) is not guaranteed: its validity depends on the 
regularization scheme adopted. If we cut all the loop-momentum components 
with a hard cutoff much smaller than the hard scale, 

|U < A uv < E x , (148) 
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The double logarithm necessarily comes from the overlap. 



34 



then the condition ( 147 ) is still valid at the loop level. As a consequence, we 
expect that the leading, double-logarithmic term of the ET shape function 
will match the QCD one. That is indeed what happens with the spatial mo- 
mentum regularization, as we have seen explicitly. On the other hand, when 
one uses a regularization such as DR or Aj_, the equality of the double-logs is 
no longer guaranteed, and indeed it does not occur in Aj_-regularization, as 
we have seen explicitly. This is because the longitudinal momentum of the 
gluon l z , or equivalently its energy e, can become arbitrarily large. For the 
latter regularizations, even for the double-logarithm, one has to come back 
to the original QCD loop diagram and perform factorization into a factor C g 



and a factor C q , as we have shown in detail. In ref. JLO] it was shown that 
the factor of 2 in the log 2 k + term in DR is a regularization effect, i.e. it can 
be removed by going to a non-minimal dimensional scheme. We explicitly 
see, with the similar A^ regularization, that by including C q the scheme- 
dependence automatically disappears. The origin of the additional factor of 
2 in the transverse-momentum regularization is related to the occurrence of 
a second double- logarithmic region for \l z \ ^> A^ (very large rapidity). 

Finally, as already noted, let us observe that in the case n 2 / we expect 
the transverse momentum cutoff to give double-logarithmic results similar 
to those from the space momentum cutoff. That is because n 2 ^ cuts the 
collinear emission at infinite rapidity. 



7 The shape function in the low-energy effec- 
tive theory 

With the A^ regularization, double logarithms are contained in C g as well as 
in C q . Since we want to confine double logarithmic effects inside the shape 
function only, let us consider from now on the As regularization only. The 
factor C q is short-distance dominated in the latter regularization, so it is 
computed once and for all in perturbation theory and "leaves the game" . 

Let us therefore return to formula ( |1 10|) for C g . Calling e = | Z |, and 
t = 9 2 , C g can be written as 
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where we have assumed As < O (Ex) and we have used the approximation 
1/ (2k + + et) ~ 8 (et - 2k + ) / (et), which is valid in DLA. This form helps 
visualizing the origin of the double logarithm. We see that contributions 
come from soft regions, where e ~ 0(k + ), as well as from hard regions, 
where e ~ O (As). In order to separate them, the simplest way is to introduce 
another UV cutoff Aet, this time well below the hadronic energy Ex (the 
hard scale of the process), such as 

k + < Aet < A s . (150) 

We can write 

C g (A s , k+) = 5Z (A s , A ET , K) + 5F ET (A ET , K) , (151) 



( )Z(;\,.;W.M = - / AS - t ld ±e(et-k, 



where 



A ET e io f 
1 1 2 As 1 , 2 A-ET ,-, p. N 

= "2 l0g k~ + + 2 l0g IT (152) 

is a coefficient function and 5F ET (A E t, k+) is the one-loop contribution 
to the light-cone function 5F ET , multiplied by the propagator: 5F ET = 
fipET i (_^ + _|_ ^o), as defined in eq. 



SF ET (A ET , k+) = - f 

Jo 



dt r 1 


dt 


e Jo 


T 


(A ET 




V k + 


)• 



- >e 2 (^J- < 153 ' 

Note that 5F ET depends only on the two scales k + and Aet- This is in 
line with the idea of a simple low-energy effective theory, which describes 
infrared phenomena characterized by the scale k + , apart from the UV cutoff 
that enters through loop effects. 
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We assume that long-distance effects can be traced by the growth of the 
coupling constant in the proximity of the Landau pole, and that the coupling 
constant must be evaluated at the transverse momentum squared [^(J: 

a s -> a s (k\) , (154) 

where 

k\ S eH. (155) 

From the expression of 5 Z we see that transverse momenta have a lower 
bound given by 

l± > l±mia = A ET k+. (156) 

According to our criteria, non-perturbative effects are absent from Z as long 

as 

^_i_min ^ Aqcd- (157) 

According to eq. (|156|) , this occurs when k + is non- vanishing, as it is for 
example if 

k + ~ O (Aqcd) , (158) 

as expected from Fermi motion (since Ket 3> Aqcd)- However, by taking 
the imaginary part of T^ u to obtain W^, i.e. the rate, the product of factors 
is converted into a convolution over k + and the point k + = is included 
in the integration range. This implies that transverse momenta down to 
zero contribute to the coefficient function in , i.e. that factorization of 
short- and long-distance effects breaks down at this point. The breakdown 
is related to the fact that we are cutting the energies of the gluons, but 
not the emission angles, which can go down to zero, implying the vanishing 
of the transverse momenta. That is one of the most important outcomes 
of our analysis. However, we believe that these long-distance contributions 
are suppressed. Let us present a qualitative argument. As we can see from 
inequalities fll56j ) and ( |157|) , transverse momenta of the order of the hadronic 
scale occur in Z for a very small slice of values of k + , 

K < ^ « Aqcd- (159) 

If the integrand is not singular in this small slice, as it is natural to assume, 
it gives a reasonally small fraction of the total. Note that the usual factoriza- 
tion of mass singularities is instead "exact". If we consider for example the 
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moments of DIS cross-section, factorization involves a splitting of the long- 
and short- distance contributions of the form 

M N (Q 2 ) = f dx B x^- 1 a DIS (x B ,Q 2 ) (160) 
Jo 

f Q2 dll A f Q2 dl 2 ± \(, f A2 dll 

= 1 + <y N a s / = 1+ 7jv a s / 1 + j N «s / 

J m 2 tj_ y J A 2 L _L J \ Jm? L J_ 

where m is the mass of a light quark. 

After the last step ( |152j ), the forward hadronic tensor takes the final form 

T QCD = ^F(k + f CD (161) 

1 — [1 + a C h \ [1 + a C q ] [1 + a 6Z] [l + a 8F ET ] , 



2v ■ Q -k+ + i0 

where the various factors have been introduced in eqs. (|77|), (ff8|), (|87|) and 
( |102|) . Taking the imaginary (absorptive) part, according to the optical the- 
orem (|44l), we have for the multiple convolution 

Wau = — / dki dk 2 dk% dk 4 5 (k + — ki — k 2 — k 3 — k 4 ) 
2v ■ Q J 

[5 (k x ) + ac h (ki)) [5 (k 2 ) + ac q (k 2 )] 

[5 (h) + adz (k 3 )} [5 (h) + adf ET (k 4 )] , (162) 

where 

f ET (k + , A ET ) = 6 (k + ) + a 5f ET (k + ) + 0(a 2 ) (163) 

is the shape function, defined in eq. (|TUD, for an on-shell quark (k' + = 0); 
moreover, we have defined 



c h (k + ) = --Im 

7T 



' C h (k + - zO) 



-k+ + iO 



= 6{k + )C h (k+)-^(-^jlmC h {k + -i0) (164) 

and analogously for the other factors^]. Typically, by taking the imaginary 

28 In ref. @, formula (9) should be replaced by f QCD {k + ) = / dki dk 2 5(k+ - fci - 
k 2 ) {8(h) + aSz(k 1 )) f ET (k 2 ), where Z= 1 + aSZ. 
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parts, for double-logarithmic contributions, we have 



5(k + )\og 2 (-k + ) + 26(-k + ) 



7T -/.• . • /'() ' • 7.\ 

d 



d/c-i 

and for single-logarithmic ones 



(-0 (-A^.) tog 2 (-*+)) (165) 



1 log (k + - iO) . 1 

s (k+) \og(-k+) + e (-k+) 



7T // . - /() ' ••/.-. 

= ^-(-9(-k + )]og{-k + )). (166) 

The last members of the above equations have to be interpreted as distribu- 
tions. In DLA, according to eq. ( |153|) , f ET up to one loop reads 

f (k + ,A ET ) = 5 (k+) + a 6 (-k+) -8(k + )\og 



M* + ) + ^("(-* + )log 2 (^)). (167) 



"A^_[_ 2 V /c_i 

2 / A £3 

2d/.+ V ' '~° g V-Jfe-. 
7.1 Evolution 

Taking a derivative with respect to the logarithm of the cutoff, we obtain 

df{k + ,A ET ) fA ET \ 6(-k + ) 

—— = - a 5 {k + ) log — — + a 

d\ogA ET \-*+/ ~ k + 

= J_(, ( _ Mlog (^i)). (168) 

Comparing the above equation with the evolution equation for the shape 
function 

d {f + g ^ ] = j dK K S (k + - k' + ; A ET ) f (k' + , A ET ) , (169) 

and taking into account that, at lowest order in as, / (k' + , A E t) = 5 (k' + ) 
holds, we find for the evolution kernel at one loop 

Aet \ e(k' + -k + ) 



K s (k + - k' + ; A ET ) = -a 5 (k' + - k+) log (^^) + a 



39 



d 

a- 



dk A 



(170) 

i - fc +) 5(k > k) [ Aet d (!± - '+) 

If we consider the Aj_-regularization, the evolution kernel for the shape func- 
tion is instead given by (eq. (|140|) ): 

/ A i \ 6(k'-k + ) 
K ± (k + -k' + ;A ± ) = -2a 5 (k' + - k + ) log ) + ^ ^ 

- 2a ik{ e{k '+- k+ MA))- (171) 

We notice that there is a factor of 2 between the kernels ( |170| ) and (|171|) for 
the shape function in the two regularizations. The kernel in DR is the same 
as that in eq. ( |171| ), with Aj_ — > \i. 

There is a clear analogy of the evolution of the shape function with the 
Altarelli-Parisi evolution equation, but with an important difference: the 
evolution kernel in this case explicitly depends on the cutoff Aet of the 
bare theory or on the renormalization point [i if we consider the renormal- 
ized theory^. All this is related to the fact that the Altarelli-Parisi evolu- 
tion involves a single collinear logarithm for each loop, while our problem 
is double-logarithmic. Let us discuss this point with a simple analogy. The 
Altarelli-Parisi evolution, or in general the usual Callan-Symanzik evolution, 
is analogous to a first-order differential equation, which is autonomous (i.e. 
time-independent) : 

! = *<*), (172) 

or, in discrete form, 

x n+1 = O (x n ) , (173) 
where O is a generic operator, such that the formal solution reads 

x n = O n (x ) . (174) 



2') 



We thank S. Catani for a discussion on this point. 
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The evolution in eq. ( |169| ) is instead analogous to an evolution equation of 
the form 

^ = h(x,t), (175) 

or, in discrete form 

x n+1 = O n (x n ) . (176) 

In the latter case there is a different evolution operator at each stepQ. 

We clarify at this point a discrepancy of a factor of 2 in the evolution 
kernel K of the shape function, computed at one loop in DR in both refs. 
|| and ||. We agree with ref. ||, where the kernel is computed from the 
Green function in the ET taking a \i derivative, as in eq. ( |171| ). We disagree 
with ref. 0, where the kernel is computed by taking the difference of the 
QCD Green function with the ET Green function and then differentiating 
with respect to //; their kernel is two times smaller than the one in eq. ( |171| ). 
The latter authors give for the QCD amplitude, in our notation, the result 

^(M 0CT = ^W-iWf^Y (178) 



-k+ + i0 V 2/ ° \k + -iO, 
They find a dependence on the renormalization point /x, which we do not 



find as the QCD diagram is ultraviolet - as well as infrared - finite |10[ . If 
we replace in their renormalization condition, which determines the kernel, 
our /i-independent result for the QCD Green function, 

F ^ QCD -^ {-\)^{>^)' (179) 

we find a vanishing kernel 0. Since the effective theory is UV-divergent 
and consequently /i-dependent, we believe that there may be a problem with 

30 In double-logarithmic problems, one can obtain an autonomous differential equation 
at the price of having a second-order equation, i.e. of the form 

5 = MX). (177) 

This, anyway, is not an evolution equation. 
31 In eq. (179) we have assumed Ex ~ O (mt). 
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the renormalization conditions. Schematically, the matrix element of a bare 
operator is of the form 



(p\0 B \p) = 1 + c-2— (p 2 y e + (finite for e -> 0), (180) 

where c is a numerical constant, p 2 refers to an overall momentum scale in 
the external state, and (p 2 ) e comes from the one-loop integral in D = 4 — 2e 
dimensions; a B im is the bare coupling of the original D-dimensional theory: 
for D < 4 it has a positive mass dimension 4 — D = 2e , and it must be kept 
fixed as we vary /i, which is just an arbitrary mass scale: 

-^-a d B im = 0. (181) 

This implies the well-known condition 

A-(p\O B \p) = 0. (182) 

One usually introduces an adimensionalized bare coupling multiplying a d ^ m 
by fi 2e , where \i is just an arbitrary mass scale as we said before, 

af m = /i" 2e a| im , (183) 

so that the bare Green function reads 

adim / 2\ e 

(p\0 B \p) = l + c^—(^ \ + (finite for e - 0) 

adim 2 

= l + c^ + ca B dim \og^ + --- (184) 
e p l 

In the minimal-dimensional scheme (MS), we include the pole term in the 
renormalization constant 



Q,adim 

Z MS = l + c^—, (185) 



and the remaining terms in the matrix element of the renormalized operator, 
(p\0 M s\p) = 1 + caf log ^ + " " " , (186) 
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since Ob = Z Or. It is only after this splitting that a dependence on fj, is in- 
troduced separately in the renormalization constant and in the renormalized 
operator^. 

The anomalous dimension is computed from the renormalization constant 
keeping a^ m fixed: 

^d^ = d(<^\ = _ 2ca ^ (187) 



d log n d log fi 

It seems to us that a vanishing kernel or anomalous dimension in the effective 
theory is obtained in ref. || because the renormalization constant Z has been 



identified with the whole matrix element (|184|) 



8 Conclusions 

We have discussed the properties of decays of heavy flavour hadrons into 
inclusive hadron states X with an invariant mass nix small compared with 
the energy Ex, rnx <C Ex- An explicit factorization procedure has been 
introduced, which holds on a integral-by-integral basis. It is based on: 

• the Cauchy theorem: it is exact and leads to the replacement of ordi- 
nary propagators with eikonal propagators in loop integrals; 

• the lowering of a hard UV cutoff from A uv > nib to A uv = A ET ^ 
nib, where Aet is the UV cutoff of the low-energy effective theory inside 
which the shape function is defined. 

This technique has led us to a clean separation of all the perturbative 
and non-perturbative contributions. We have found that, while the exact 
kinematics of the original QCD processes involves a Wilson line off the light- 
cone for the final light quark, in the low-energy effective theory the light 
quark is necessarily described by a Wilson line on the light-cone. 

We have analyzed the shape function / (k + ) to find out which long- 
distance, non-perturbative, effects are contained in and which are not, in 
different regularization schemes. We found that f(k + ), contrary to naive 

32 In the notation of ref. §, log(O s ) = dbxf B (f) /<91og£, with £ - l/v 7 ?- 



43 



physical expectations, has no direct physical meaning even in double loga- 
rithmic approximation, as it represents a partial contribution to the complete 
physical process. Changing regularization, we have explicitly shown that the 
leading double-logarithmic contribution to / (k+) can be changed by a fac- 
tor of 2, i.e. that the shape function is substantially regularization-scheme 
dependent. Only after summing the shape function with the other contri- 
butions, is a physical, scheme-independent result recovered. We have also 
shown that in lattice-like regularization the shape function factorizes a large 
part of the non-perturbative effects: it contains all the double logarithmic 
contributions of the full QCD process. 

Subtracting from the forward hadronic tensor T® v , step by step, each 
of the perturbative components, we end up with an explicit representation 
of the perturbative and non-perturbative effects. For instance, at one-loop 
order in DLA we have the result (see eq. (TE2")): 



(188) 

The coefficient Ch is a hard factor that takes into account the fluctuations 
with energy e in the range Ex < e < ttib- The other two coefficients, C q 
and SZ, are short-distance-dominated in lattice-like regularization schemes; 
5F ET is long-distance-dominated in any regularization. The tensor W^ u , i.e. 
the rate, is obtained (as usual) by taking the imaginary part of T^ u . 

Another outcome of our analysis is that, contrary to single logarithmic 
problems, factorization in this (double- logarithmic) problem is not exact, 
even in lattice-like regularization schemes. Some long-distance effects are 
present in the coefficient function: they come from gluons with a large energy 
but with a very small emission angle and consequently with a small trans- 
verse momentum. These non-perturbative effects in the coefficient function 
however are expected to be suppressed on physical grounds, as they occur in 
a small region of the phase space for a moderately large cutoff of the effective 
theory. 

Finally, we have clarified some discrepancy in the literature about the 
evolution kernel for the shape-function computed in double logarithmic ap- 
proximation inside dimensional regularization. 
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Figure 1: Pictorial description of factorization in the effective theory. 
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Figure 2: Vertex corrections to the light-cone function F® (k + ) 
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c c s c h 

Figure 3: First decomposition in T^ u . The thick, thin and double lines 
represent the massive quadratic, massless quadratic and time-like eikonal 
propagators, respectively. 
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€ — plane 




Figure 4: Poles of C s in the Zo-plane. The crosses labelled by b, u and g 
represent the particle poles in the static beauty, up and gluon propagators, 
respectively. The crosses labelled by u and g, instead, represent the antipar- 
ticle poles in the lip-quark and gluon propagators. 
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Figure 5: Second decomposition in T^ v . The thin and double lines represent 
the massless quadratic and eikonal propagators, respectively. The symbols 
(v) and (n) indicate that the eikonal propagators must be taken along the 
directions v and n. 
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